Polaron relaxation in self-assembled quantum dots: Breakdown of the semi-classical model 
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We calculate the lifetime of conduction band excited states in self- assembled quantum dots by taking into 
account LO-phonon-electron interaction and various anharmonic phonon couplings. We show that polaron 
relaxation cannot be accurately described by a semi-classical model. The contributions of different anharmonic 
decay channels are shown to depend strongly on the polaron energy. We calculate the energy dependence of 
polaron lifetime and compare our results to available experimental measurements of polaron decay time in 
InAs/GaAs quantum dots. 
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Carrier relaxation in quantum dot (QD) is of fundamen- 
tal importance for potential applications {e.g. laser, photo- 
detectors...). In self-assembled quantum dots, the energy 
separation between the ground and first excited states is of 
the order of tens of meV. Early theoretical studies of car- 
rier relaxation in these nanostructures led to the so-called 
phonon bottleneck, in which the phonon-assisted relaxation 
is strongly inhibited because of the discreteness of electronic 
states Further studies predicted the existence of a strong 
coupling regime between electron and longitudinal optical 
(LO) phonons in self-assembled QDs^. This strong cou- 
pling regime leads to the formation of polarons, which are 
entangled electron-LO-phonon quasi-particles. Polaron states 
were first observed in intra-band magneto-transmission exper- 
iments done on n-doped QDs^'^'^. More recently, evidence of 
a strong coupling regime has also been obtained for holes^ and 
electron-hole pairs^ confined in QDs. A striking consequence 
of the polaron formation is that electrons can no longer re- 
lax by LO-phonon emission. Nevertheless, efficient polaron 
relaxation (lifetime Tpoi of few tens of picosecond) has been 
measured in QDs doped with one electronic and more re- 
cently with two electronsi^. 

So far, there are two models in the literature to describe the 
polaron relaxation in QDs. Li et al^^ were the first to con- 
sider the carrier relaxation as triggered by the intrinsic (bulk 
related) instability of LO-phonons. They discussed a semi- 
classical model in which a phonon damping was phenomeno- 
logically added when solving the time-dependent Schrodinger 
equation for the coupled electron-LO phonon system. They 
were able to derive a formula for the polaron lifetime, which 
depends on two parameters: the strength of the electron-LO- 
phonon coupling and the damping rate. The later were taken 
as fit parameters in Refs 10 and 11 to successfully interpret 
the increase of Tpoi with the polaron energy Epoi in the 40-52 
meV interval. However, this model fails in explaining the de- 
crease of Tpoi for polarons with higher energy (up to ~ 60 
meV). On the other hand, Verzelen et aM and Jacak et aM 
have applied the Fermi golden rule to polaron states, assuming 
the anharmonic mechanism proposed by Vallee and Bogani 
for bulk LO-phonons (LO LO + TA)^^. This assumption 
leads to the existence of a narrow energy window for relax- 
ation : only polarons with an energy Epoi in the ~ 35-44 meV 
interval can disintegrate. This model also predicts an increase 
of Tpoi with Epoi > hujLO (where ujlo is the LO-phonon 



frequency), but of course fails when applied to higher energy 
polarons. 

In this paper we present a detailed study of polaron relax- 
ation in QDs. We calculate the lifetime of excited states in 
self-assembled QDs by fully taking into account their pola- 
ronic nature as well as the phonon anharmonicity. We show 
that the variation of Tpoi with Epoi is not governed by the sole 
weight of its LO-phonon component (like in the work of Li 
et al), neither can be restricted to only one bulk disintegration 
channel (like in the works of Verzelen et al and Jacak et al), 
but that different anharmonic channels have to be taken into 
account depending on Epoi . As we show below, the polaron 
lifetime is equal to the product of its LO-phonon weight by 
the decay rate of a bulk LO-phonon which would have the po- 
laron energy. We obtain a good overall agreement with the 
experimental data for any polaron energy. In particular, the 
model allows to explain the measured unexpected decrease of 
Tpoi with increasing Epoi at high energies. 

We consider a quantum dot hosting one electron in its con- 
duction band (as obtained e.g. by doping). The system is 
modelled by the following Hamiltonian: 
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where He is the purely electronic Hamiltonian, Hyib the vi- 
brational one and He-ph represents the different electron- 
phonon interactions. The vibrational Hamiltonian (assumed 
to be bulk-like), can be expressed as Hyi^ = Hph + K 
where Hph is the harmonic hamiltonian of phonons while Va 
is the anharmonic part. In actual self-assembled QDs, the en- 
ergy spacing between ground and first excited level is much 
larger than the energy of acoustic phonons which have the QD 
size wavelength. Therefore, we can neglect direct couplings 
related to acoustic phonons and consider only the Frohlich 
Hamiltonian: 



E 



(2) 



where Cf 



ulo and uto are 



respectively the zone-center LO and TO modes frequencies, 
£00 is the high frequency permittivity and Vcr is the crystal 
volume. 

Let us define Hq = He -\- Hph + Hf. The polaron states 
are eigensolutions of Hq. We recall briefly in the following 
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FIG. 1: (a) Full lines: Polaron energies as a function of the electronic en- 
ergy separation between the \s) and \p) levels. Dashed lines: energies of 
non-interacting |s, 1} and |p, 0} level, (b) Scheme of the interaction between 
electronic levels and LO-phonon, leading to the polaron formation. Anhar- 
monic coupling with the phonon reservoir and the far-infrared absorption are 
indicated (see text). 



their mainaspects (for more details on polaron eigenstates see 
e.g. Ref 6). The polaron wavefunctions are entanglements 
of decoupled states |a, n) where a = s^px or py is the elec- 
tron wavefunction (solution of H^) and n is the number of 
LO phonon modes (solutions of Hph) entering in the forma- 
tion of the polaron. Since the LO-phonons that are involved 
in the Frohlich couplings have small wavevectors^, we can as- 
sume that LO phonon are dispersionless. Within this approx- 
imation, the Frohlich interaction between the discrete zero 
phonon level |p, 0) and the one phonon flat continuum |s, Iq) 
can be treated as an interaction between two discrete levels, 
namely |p, 0) and |s, Igp), where the phonon mode \lsp) is 
defined by: 
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where \lq) are bulk LO-phonon states, and 



sp 



1(5, 0)p is the coupling strength between 

|p, 0) and |s, Isp). The two resulting polaron states (i=l,2) 



can be written: 



(4) 



The polaron energies Ei are solutions of {Ei — hui^o){Ei — 
Esp) = Vgp, where Egp is the energy separation between the 
the \s) and \p) levels (experimentally, different Esp values 
correspond to QDs with different sizes and/or compositions). 
The energy origin is taken at the |5, 0) leveP^, so that the po- 
laron energy Ei corresponds to the one of the far-infrared pho- 
ton (see Fig 1). The weight of the one-phonon components of 
the polaron states i can be expressed as a function of the po- 
laron energy Ei as: 



{Ei - hco] 
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The anharmonic perturbation Va triggers the polaron relax- 
ation. It couples the polaron state via its LO-phonon compo- 
nent to a reservoir of multi-phonon states. Using Fermi golden 
rule, the broadenings of the polaron states are given by: 
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FIG. 2: (Color online) Calculated F^^ versus the energy: Full calculation 
(solid line); LA-i-LA channel (Klemens channel) in dashed line; LA-i-TA chan- 
nel in dotted line; LA-i-TO channel in dashed-dotted line; TA-i-LO channel in 
dashed-double-dotted line (T=OK). Inset: Schematic of the different anhar- 
monic processes. 



ri(Ei) = lftl'rp''(^i) (6) 



V^\E,) = 2T:Y^\{v\Va\\sp)\^KEi - E,) (7) 

V 

where v label the multi-phonons states into which the po- 
laron disintegrate. The \f5i\^ term follows from the fact that 
Va acts only on the LO-phonon part of the polaron eigen- 
state. Note that the semi-classical approach leads to Vi{Ei) — 
lAprP^(ncjLo), where n/rP^(ncjLo) is the bulk phonon life- 
time. In our model instead, the polaron lifetimes are evalu- 
ated at the polaron energies. Since the optically probed state 
has energy Ei that can significantly differ from fic^Lo, the an- 
harmonic coupling strength as well as the density of multi- 
phonons final states \v) resonant with Ei may greatly differ 
from the bulk. Therefore, as we show below, rP'^(£^^) can be 
very different from V^'^ihwi^o). 

Linewidth of bulk LO-phonon r(ficjLo) has been calculated 
in several works: first principles calculations have been used 
by Debernardi et al}^'^^, while Klemens^^ or more recently 
Barman and Srivastava^^ have used a more phenomenological 
approach. Here, we generalize the later model to study the 
relaxation of polarons in QDs. We consider only the cubic 
part of the anharmonic vibrational hamiltonian^^: 
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where ft is the unit cell volume, Ncr = Vcr/^ and 
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Va describes the coupling of phonon of branche jo with 
two-phonon states of branches ji and j2- Each (ji,j2) cor- 
responds to a particular disintegration channel. The 3 ! factor 
accounts for equivalent processes obtained by cyclic permuta- 
tion of the labels in the triple summation. 

As the \lsp) LO-phonon mode is formed of bulk phonons 
with small wavevectors, momentum conservation in Eq.[8]im- 
plies a decay into a pair of phonons with opposite wavevec- 
tors. Therefore, one gets approximatively at low temperature: 
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(10) 



Following reference 21, which extents Klemens' workiS, 
we use: 

|P(/co, jo;^i, ji;^2, j2)| = ^^fco,jo^fci,ii^fc2,j2 (11) 

where 7 is the mode-averaged Gruneisen's constant (we use 
7 = 1.3 in the following^^'^^), p is the mass density and c is 
the average acoustic phonon speed (3/c = 1/cla + 2/cta)- 
T^^{E) can be decomposed into the contributions due to the 
different channels (ji , ) • 



(12) 



Using Eq.[TTl we can express each of these contributions as: 
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As in Ref 20, we model acoustic modes with Debye's 
isotropic continuum scheme. The linewidth induced by a 
given channel (ji, ^2) can be expressed as: 



" (J1,J2) 



(14) 



where Cj is the sound velocity of branch j for acoustic 
branches {cj^ = duj^/dkj^ for optical branches). In Eq[T4l 
ujj-^ and cjj2 are the phonon frequencies that satisfy both deltas 
in Eq[T3l For a non-zero temperature, this linewidth has to be 
multiplied by (1 + J (1 + rij^), where Uj. is the Bose occu- 
pation factor at the energy hjj. . Figure[2l shows the calculated 
PP^ versus the polaron energy, as well as the contributions to 
pph from different channels: F^^^ ^j^y F^^^ ^^^^ (Klemens 

channel), T^^^ (Vallee-Bogani channel) and E^^^^^^^. 

We obtain very strong variations of FP'^ within the range of 
polaron energies that can be obtained in InAs/GaAs QDs. 
This radically differs from the constant value FP^(ficjLo) of 
the semi-classical model. In order to understand the origin of 
these variations, let us consider the calculation of the Klemens 
channel : 
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FIG. 3: (Color online) Polaron decay time as a function of the polaron en- 
ergy: measurements from Ref 11 in symbols, full calculation (presented in the 
text) in solid line, calculation involving only the Klemens channel (LA+LA) 
in dotted line and best fit obtained with the semi-classical model in dashed 
line (Vsp(fit) = 8meV, tlo = 15ps). 



(LA,LA) 



(15) 



The dependence arises from two contributions: the den- 
sity of final states with two LA phonons with energies equal to 
half of the polaron one increases like E'^ (in the Debye model 
we are considering), while the anharmonic coupling strength 
also increases like E'^. The disintegration rate F^^^ 
increases with E until the energy of the emitted LA phonons 
reach their maximum value (corresponding to the emission of 
two zone-edge LA phonons, for E b2 meV). We discuss 
further the implications of these results below. Let us first 
point out that we have checked that the Fermi golden rule can 
be safely applied in our case: (i) the anharmonic self-energy^^- 
(the Hilbert transform of ^^^{E)) is negligible (on the 0.1 
meV energy scale, except at channel edges) with respect to the 
typical variation of E studied here, (ii) dV^^/dE <C 1 except 
at channel edges, so that the time decay of the polaron popu- 
lation is exponential. Finally, as shown in Fig 2, we find that 
the (LA,LA) channel is the dominant mechanism in bulk^^. 

The final polaron lifetime Tpoi = ^/E^ (including both the 
results in Fig 2 and the energy-dependent weights |/3ip) as a 
function of polaron energy Epoi = ^i=i,2 is shown in Fig 3 
(solid lines), and compared to the data of Ref ITTI (symbols). 
We observe a good overall agreement with the measurements 
for any polaron energy. Before discussing the complex en- 
ergy dependence of Tpoi, let us consider the applicability of 
the semi-classical model. Firstly, the best fit using this model 
(dashed line in Fig 3) is obtained for Vspi^fn) — 8 meV, which 
is about twice the calculated value of the Frohlich coupling 
strength for the QDs studied in Refs JO, and JJj. Note that 
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Vsp ^ 4 meV is found to be roughly independent of the po- 
laron energy, so that the factor of two in Vsp(fit) cannot re- 
sult from the inhomogeneous distribution of QDs. Using in- 
stead the calculated value of Vsp in the semi-classical for- 
mula would lead to a much steeper increase in Fig 3, since 
Tpoi ^ h{Epoi - hu;Lo)^/[^^^{^Lo)V^p] for high energy 
polarons^^. Secondly, the semi-classical model predicts a con- 
tinuous increase of Tpoi with the absolute value of the energy 
detuning \Epoi — hwLo\- Our results do not support this pre- 
diction, as discussed in the following. 

Indeed, the existence of different channels, with specific en- 
ergy windows, leads to important variations of Tpoi with Epoi, 
in particular near the windows edges. We associate the pro- 
nounced discontinuity in the data of Ref [111 at :^ 52 meV with 
the high-energy limit of the Klemens channel. According to 
Fig 2, this disintegration channel dominates over a large en- 
ergy interval. This is in agreement with other data of Zibik et 
al, who measured the temperature dependence of the polaron 
lifetime at Epoi = 44 meV and found that the dominant mech- 
anism at this energy was the (LA,LA) one^^. Furthermore, 
we clearly see from Figs 2 and 3 that the dominant disinte- 
gration path depends on the polaron energy. We associate the 
important decrease of Tpoi above 52 meV with the (LA,TO) 
channel. It is worth stressing that the efficiency of the anhar- 
monic disintegration increases in this region despite the fact 
that the weight of the one-phonon component continuously 
decreases with increasing E. This counterintuitive result high- 
lights again the importance of employing an accurate descrip- 
tion of the anharmonic coupling and of the two-phonon den- 
sity of states. Finally, we predict that the disintegration of low 
energy polarons (Ei <C Hulo Fig la) should be governed 
by the (LA,TA) channel. Note that lower branch polarons 
in Fig la) with energies near (but below) Hujlo have 
lower lifetime than the upper branch ones, and also that Tpoi 
increases much faster with increasing detuning in the lower 



branch than in the upper one. Let us finally stress that we do 
not have any fit parameter in our model, but consider bulk- 
like values for the different material parameters, in particular 
an averaged channel-independent value for 7. This approx- 
imation, together with the very simple (isotropic and linear) 
phonon dispersion we consider, may explain the discrepancies 
between calculated and measured lifetimes, especially near 
the zone edges. 

In conclusion, we have presented in this work a complete 
study of polaron relaxation in self-assembled QDs. We have 
shown that the decay rate of a QD polaron is equal to the prod- 
uct of its phonon weight by the decay rate of a zone-center 
bulk LO phonon which would have the polaron energy. The 
good overall agreement of our calculations with reported data 
on the population lifetime of excited QD states strongly sup- 
ports the mandatory use of an accurate description of phonon 
anharmonicity within the polaron framework to describe the 
relaxation process in these structures. Moreover, we have pre- 
sented a critical discussion of the previous models proposed 
to tackle the polaron relaxation, in particular the often quoted 
semi-classical one. Finally, we hope that our work will stimu- 
late further experiments, in particular for low-energy polarons 
pertaining to the lower branch (see Fig la), which are ex- 
pected to be less sensitive to the anharmonic disintegration. 
It would also be interesting to study higher energy states, for 
which higher order anharmonic terms may become important. 
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